Quantum gates and multipartite entanglement resonances realized by motion 
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We demonstrate the presence of genuine multipartite entanglement between the modes of quantum 
fields in non-uniformly moving cavities. The transformations generated by the cavity motion can be 
considered as multipartite quantum gates. We present two setups for which multi-mode entangle- 
ment can be generated for bosons and fermions. As a highlight we show that the genuine bosonic 
multipartite correlations can be resonantly enhanced. Our results provide fundamental insights into 
the structure of Bogoliubov transformations and suggest strong links between quantum information, 
quantum fields in curved spacetimes and gravitational analogs by way of the equivalence principle. 
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I. INTRODUCTION 

During the past decade research in the area of rel- 
ativistic quantum information has addressed questions 
concerning the inherently relativistic aspects of quantum 
phenomena, unveiling the close connections between ef- 
fects in quantum field theories and quantum information 
theory. 

In this context several models to store and access bi- 
partite quantum correlations in relativistic settings have 
been proposed, including Unruh-DeWitt type detector 
models [1], free field modes [2-4] and cavity modes (see 
Refs. [5-10] for details and Ref. [11] for a recent review). 

In these situations bipartite entanglement was exten- 
sively studied but, save for a few exceptions [12-16], the 
role of genuine multipartite correlations in these studies 
was marginal. However, not only are multipartite corre- 
lations expected to feature in relativistic scenarios, e.g., 
in the transformation of Minkowski modes to Rindlcr 
modes [2], their possible realizations even grow vastly 
with increasing size and complexity of the systems. More- 
over, the classification of genuine multipartite entangle- 
ment will allow for more detailed characterization of rel- 
ativistic effects and may be used for high-sensitivity tests 
of the "quantumness" of correlations. In fact, the identi- 
fication of quantum correlations can be a keystone to the 
experimental verification of effects in quantum field the- 
ory, such as the Hawking effect in analog fluid systems, 
see Ref. [17]. 

Here we want to focus on the generation of genuine 
multipartite entanglement by the non-uniform motion of 
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cavities that contain relativistic quantum fields. We em- 
ploy the techniques developed in Refs. [G, 7], for bosonic 
and fermionic fields respectively, where the moving cavi- 
ties follow worldtubes that are composed of segments of 
inertial motion and uniform acceleration. In this frame- 
work we ask the questions: Can the Bogoliubov trans- 
formations that are induced by the non-uniform motion 
generate genuine multipartite, quantum correlations? If 
this is the case, can we quantify and/ or classify the aris- 
ing entanglement? 

We present two scenarios for which we can indeed cre- 
ate, and partially classify, such correlations, thereby real- 
izing quantum gates by motion. First, in Sec. II, we con- 
sider three individual cavities, labelled Alice, Rob and 
Charlie, that share pairwise bipartite entanglement in 
their initially common rest-frame, before Rob's cavity un- 
dergoes non- uniform motion, see Fig. 1. In the second 
scenario we consider Rob's cavity on its own in Sec. III. 
In these scenarios the significantly different advantages of 
bosonic and fermionic systems become apparent. In the 
bosonic case the transformations induced by the cavity 
motion generate multi-mode entangling quantum gates. 
Moreover, we show that the genuine multipartite char- 
acter of the bosonic entanglement can be enhanced res- 
onantly by appropriate timing of the cavity's trajectory 
segments. The qubit structure of the fermionic systems, 
on the other hand, allows for a clear classification of the 
arising multipartite correlations. The cavity motion ef- 
fectively acts as a multipartite quantum gate producing 
Dickc states and W states. We present explicit pure state 
decompositions for these classifications in Sees. II B and 
IIIB, respectively. 

The calculation of the Bogoliubov transformations for 
the non-uniform cavity motion are performed pcrturba- 
tively in terms of a parameter h that physically repre- 
sents the product of the (rest-frame) cavity width S and 
the acceleration at the centre of the cavity. Using the 
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Bogoliubov coefficients obtained in this way we can per- 
form all other computations analytically. For simplic- 
ity we restrict our considerations to (1 + l)-dimensional 
Minkowski spacetime where the metric tensor rj^ has the 
signature (—,+)• The setup can be naturally extended 
to accommodate additional spatial dimensions that enter 
into the effective mass of the field via the correspond- 
ing transverse momenta [6] . We use units where Planck's 
constant and the speed of light are dimensionless, i.e., 
h = c = 1. By 0(x) we denote a quantity for which 
0{x)/x remains bounded as x goes to zero. We use an 
asterisk and a dagger to denote complex and Hcrmitian 
conjugation respectively. 



II. MULTI-CAVITY ENTANGLEMENT 

As the first scenario we consider three cavities, Alice, 
Rob and Charlie. Each cavity represents an individ- 
ual spacetime in which a rclativistic quantum field is 
confined by appropriate boundary conditions. We can 
assume without loss of generality that all three cavities 
arc manufactured in the same way and that they are 
initially (in the "in-region") at rest with respect to 
each other. At t = Rob's cavity starts to accelerate 
linearly and uniformly. We assume that the worldtubc 
of Rob's cavity consists of periods of inertial motion and 
uniform acceleration such that the cavity remains rigid 
throughout the journey, i.e., the cavity length <5 > 0, as 
measured by the co-moving observer, is constant. For 
the perturbative treatment we further assume that the 
accelerations in all individual segments are small with 
respect to the inverse cavity length 1/5. Fig. 1 shows 
the spacetime diagram of this setup for a sample travel 
scenario. 



A. Scalar fields 

Let us consider a real, scalar field <\> of mass m > 
in Rob's initially inertial cavity in the "in-region" . The 
field satisfies the Klein-Gordon equation (— 0+m 2 )(f) = 0, 
where □ is the scalar D'Alambertian. We follow the pro- 
cedure laid out in Ref. [6] and represent the confinement 
to the cavity by imposing the Dirichlct boundary condi- 
tions <p(t < 0,a) = 4>(t < 0, 6) = 0. We thus obtain a 
discrete spectrum of mode functions (f> n in the in-region, 
where the field can be decomposed as 

<j> = ^ (<t> n a„ + 4>* n 4 ) • (!) 
n 

The annihilation and creation operators, a n and a* n , sat- 
isfy the usual commutation relations [a„,aj„] = S nm , 
The vacuum state of the corresponding Fock space is an- 
nihilated by all a n 's, i.e., a n \ 0) =0. 
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FIG. 1. Space-time diagram of the multi-cavity setup: While 
Alice's and Charlie's cavities remain at rest Rob's cavity un- 
dergoes non-uniform motion. Rob's cavity of width 8 — (b — a) 
follows a worldtube that consists of alternating segments of in- 
ertial motion (red) and uniform acceleration (green). The in- 
dividual segments in this example travel scenario are of equal 
duration t/2 and equal acceleration h/S as measured at the 
centre of Rob's cavity. The dashed lines indicate the light 
cone at (t, x) = (0, 0). 



Before we construct our initial state let us briefly recall 
the characteristics of genuine multipartite entanglement. 
Any pure quantum state that can be written as a ten- 
sor product with respect to any bipartition is called bi- 
separable. Generalizing this notion to mixed states, any 
mixed quantum state is called bi-scparablc if it admits 
at least one decomposition into a convex sum of pure, 
bi-separable states. Conversely, all states that are not 
bi-separable are called genuinely multipartite entangled. 
See, e.g., Refs. [18, 19] for more details. 

We now proceed by constructing an initial state that 
contains no such genuine multipartite entanglement. We 
select two of Rob's modes, k and k' , and entangle them 
with modes A and C in Alice's and Charlie's cavities 
respectively. The initial, bi-separable state of these four 
modes we denote by p%j C = |$ ± ) ($ ± |, with \ <& ± ) = 
I ) Ak I ^ ± )h'Ci wnere the bipartite entangled states are 
given by 

l^ ± U = 7 5 (|0>±|lA>|l fc >), (2a) 
l<^W = 7f(|0}±|l fe '}|lc)). (2b) 

At t = we start to accelerate Rob's cavity uniformly and 
linearly and we let it follow a worldtube that consists of 
segments of inertial motion and uniform acceleration. An 
example trajectory is shown in Fig. 1. After an arbitrary 
number of such trajectory segments we assume without 
loss of generality that the cavity remains inertial in the 
"out-region". The mode functions <j) m and their corre- 
sponding annihilation and creation operators a m and a\ n 
in the out-region are related to their in-region counter- 
parts by a Bogoliubov transformation with coefficients 
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(%mn and /3 mni i.e., 

4>m = y^(amn0n + flmnfin) ( 3 ) 
n 

and 

= X] ( a ™« a « _ ■ ( 4 ) 

n 

In the case where the length S = (b — a) of the rigid cavity 
is fixed and the individual accelerations of the cavity at 
any point of the trajectory are small compared to 1/5 the 
Bogoliubov coefficients can be computed analytically as 
a Maclaurin expansion, i.e., a mn = a^n + «mn + 0{h 2 ) 
and p mn = + pSb + 0(h 3 ), see Ref. [6]. Here the 
superscripts indicate the power of the expansion parame- 
ters hi , which represent the products of the cavity width 
and the acceleration at the center of the cavity in the i-th 
trajectory segment. For ease of notation we will suppress 
the subscript i in the remaining discussion. 

Using well-known standard procedures (see, e.g., 
Refs. [6] or [8]) the Fock states \ rik) and \ny) of Rob's 
cavity in the initial state p^n C can be transformed to the 
out-region basis | hi ) to obtain the transformed density 
matrix p^n C - Subsequently, we trace over all of Rob's 
out-region modes except k and k' , gaining the reduced 
density operator p^kk'C of the four bosonic modes A, k, k' 
and C . Keeping terms up to first order in h we find that 
PAkk'c ^ s effectively a state of two qubits, the modes A 
and C, and two qutrits, the modes k and k'. For such 
states a general quantification of genuine multipartite en- 
tanglement proves to be cumbersome, as there is no gen- 
eral classification scheme for multipartite systems beyond 
qubits. However, we can construct an inequality [19-21] 
that acts as a witness for genuine multipartite entangle- 
ment by comparing diagonal- and off-diagonal elements 
of PAkk'C- Since the diagonal elements that are newly 
generated by the Bogoliubov transformation are at least 
0(h 2 ) the witness inequality takes the simple form 

V2\ < lc | < 2 V | ( 2 fc | ( U | p% k , c I 6) I - 0(h 2 ) < 0. 

(5) 

Evaluating the matrix element we can recast (5) as 

^\&-O(h 2 )<0. (6) 

Within the perturbativc regime this inequality is violated, 
i.e., genuine multipartite entanglement is detected, when- 
ever P kk \ 7^ 0. This is the case for all mode pairs of 
opposite parity, i.e., if (k — k') is odd. In fact, (6) is 
not only a witness, but also a lower bound to a mea- 
sure of genuine multipartite entanglement, the average 
minimal mutual information across any bipartition, min- 
imized over all possible decompositions, which is based 
on the Renyi-2 entropy [21]. Moreover, to first order in 
the Maclaurin expansion p^ RC is a pure state for which 
the bound is tight, as shown in Ref. [20]. In this re- 
spect the bounds provided by the witness inequality offer 



an immense advantage with respect to the direct com- 
putation (if possible) of entropy-based measures, which 
require the perturbative corrections to be calculated at 
least to second order in h, see Ref. [10]. 

Conceptually the creation of multipartite correlations 
in the scenario with three, initially pairwise, bipartite 
entangled cavities, could be interpreted in the follow- 
ing way: The non- uniform motion correlates the non- 
interacting modes k and k' in Rob's cavity, such that 
all bi-partitions of the system become entangled. As ex- 
pected, the reduced state of Alice's and Charlie's cavities 
is unaffected by the motion of Rob's cavity and, conse- 
quently, the modes A and C are still uncorrelated in the 
out-region. The reduced state of the modes k and k' on 
the other hand becomes entangled, i.e., to first order in 
h the negativity Af of p kk , = T^A,c{pAkk'c) is g iven by 
Af = |/^,|/4. The negativity 

i 

which captures how the eigenvalues Xi of the partially 
transposed density matrix fail to be positive, see Ref. [22], 
is a useful measure in the context of the perturbative cal- 
culations at hand. It allows the quantification of entan- 
glement at leading order, while entropic measures of en- 
tanglement rely on the quantification of mixedness, which 
is not altered to first order in h and thus require the per- 
turbative calculations to be extended to higher orders. 
The fact that both the genuine multipartite entangle- 
ment of PAkk'C an d the bipartite entanglement of p kk , 
are controlled by f3 kk , supports the interpretation above. 
However, we shall see that this naive view does not hold 
for the corresponding fcrmionic scenario in Sec. II B. 

Nonetheless, the entanglement of the four- mode system 
is genuinely multipartite. And, moreover, the responsi- 
ble coefficient \P kk \\ can be resonantly enhanced, within 
the limitations of the perturbative regime, by appropriate 
travel scenarios when a massless scalar field is contained 
within the cavity [10]. 

Such resonances can occur if the chosen travel scenario 
consists of TV identical building blocks of two or more dif- 
ferent, inertial or uniformly-accelerated, trajectory seg- 
ments. If the overall proper time r of one such building 
block, measured at the centre of the box, satisfies the 
necessary resonance condition 



where n is a positive integer, and the coefficient 
of a single such building block is nonzero, then the cor- 
responding coefficient after N repetitions scales linearly 
with N, i.e., 

(Mil = ^V(/3iC ■ (9) 

Sec Fig. 2 for an illustration of resonances for different 
mode pairs. This scaling is valid within the perturbativc 
regime, i.e., as long as Nh <C 1. 
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FIG. 2. The linear coefficients \p^ n \ for the bosonic Bogoli- 
ubov transformations, are shown for a real, massless scalar 
field in (1 + 1) dimensions. The travel scenario has N seg- 
ments of uniform acceleration h/5 and duration r/2 as mea- 
sured at the cavity's centre, separated by (iV — 1) segments 
of inertial coasting of the same duration (to first order in 
h), as illustrated in Fig. 1 for N = 2. The curves in 
Fig. 2 are plotted for N = 15, (m,n) = (k, k') = (1,2) 
(blue, solid) and (m,n) = (k',k") = (2,3) (red, dashed) and 
u := hr/[4:5 atanh(/i/2)]. The vertical dashed lines indicate 
the potential resonance times for (k,k') and (k',k") respec- 
tively. 



B. Dirac fields 

Let us now consider the analogous scenario for cavi- 
ties containing fcrmionic rather than bosonic quantum 
fields. In particular, let us assume that Alice's, Rob's 
and Charlie's cavities confine Dirac fields, as discussed in 
Refs. [7, 8]. In the in-region the field in Rob's cavity sat- 
isfies the Dirac equation (ij^d^ — m)ip = and the mode 
functions obey the boundary conditions 

Ml7VlM*=« = (^7 7Vc')*=& = • (10) 

Here 7 M (/x = 0, 1, 2, 3) are the usual Dirac matrices sat- 
isfying the anti-commutation relations {7^,7^} = 2r]V . 
Due to the boundary conditions the field modes can be 
labelled by an integer n and the field can be decomposed 
as 

i> = Y, b n*Pn + J2 C n^- ( U ) 

ti>0 n<0 

The operators b n and c n annihilate particles and antipar- 
ticles respectively. The vacuum state of the cavity is de- 
fined by b n I )) = Cn I )) = Vn, while the single par- 
ticle and antiparticle states, |l m )) + and ||l n )) _ , are cre- 
ated from the vacuum by bJn and c n respectively. The 
fermionic operators further satisfy {b m ,b n } = {c m , c n } = 
Smm while all other anti-commutators vanish. 

In analogy to the bosonic states (2) we select initially 
bipartite entangled fcrmionic states that correlate the 
cavities of Alice and Rob, as well as Rob and Charlie 



respectively, i.e., 

1^)^ = ^(10)) ± II U))" II U + ), (12a) 

i^»flc = 7f ± i 1 **'))" net), (12b) 

where we have chosen particle and antiparticle modes in 
agreement with charge superselection rules. Also, we se- 
lect a particular mapping from the fermionic Fock space 
of each individual cavity to the tensor product structure 
of the multi-qubit Hilbert space. We use the convention 
that the single-qubit Hilbert spaces for particle and an- 
tiparticle modes are ordered by ascending mode labels, 

c^|0))=:||l p ))-®|l g )) + ®||l,.)> + , (13) 

where p < 0, q, r > and p < q < r, such that the 
corresponding creation operators need to be ordered ac- 
cordingly. We reverse this ordering in the adjoint Hilbert 
space. For clarity we have used the tensor product sym- 
bol here explicitly, but we omit it again in the remaining 
paper. This step is necessary, since fermionic Fock spaces 
are not naturally equipped with the same tensor product 
structure as bosonic Fock spaces, see, e.g., Ref. [23] for a 
discussion of this issue. However, we have verified that 
this fermionic ambiguity does not affect the results pre- 
sented in this article. Having adopted the convention of 
Eq. 13, we continue, as previously, by transforming the 
initial state ||$ ± )) = || 4^ }} Ak || ^ ± } K , C with the corre- 
sponding density operator = |$±)) (($±1, to the 
out-region states || 0)), |l m )) + and ||ln)J~ 1° obtain the 
transformed state Qahq- The Bogoliubov transformation 
that relates the in-region modes ip n and the out-region 
modes ip m is given by 

j>m = ^ Arm^n (14) 
n 

where n G IN and the fermionic Bogoliubov coefficients 
can be expanded in a Maclaurin series in the parameter 
h as Amn = + Attn + 0(h 2 ), see Ref. [7]. 

Subsequently we trace over all particle and antiparticle 
modes except n > and k < 0, which leaves us with the 
reduced state Qakk'C- ^° nrs ^ or d er in the parameter h no 
new diagonal elements are generated in the state trans- 
formation, which makes it easy to identify an inequality 
that acts as a witness for genuine multipartite entangle- 
ment. However, due to the Pauli exclusion principle, we 
can not employ a witnesses of the type of Eq. (5) for the 
fcrmionic system. Instead we construct the witness 

r«i K ,i + ((i K i^ciiK')nic)) + i + (15) 
r((i c i + ((uii^ KK , c ii K onic)) + i - o(h*) <o, 

which we can express as 

IKk'I - o(h?) <o. (i6) 
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The witness employed in Eq. (16) is again a lower bound 
to the minimal average mutual information across any 
bipartition [21], although, this time, it is not tight for 
pure states as Eq. (5). However, there are other con- 
ceptually intriguing features appearing for the fermionic 
four-qubit state Qakk'C- First, we notice that the un- 
perturbed reduced density operators g^ c and g^, are 
both maximally mixed. This means that, in contrast to 
the bosonic case, no bipartite entanglement between the 
fermionic modes k and k! can be generated from the ini- 
tial state g^ftc by small perturbations. The negativity, 
which is nonzero for all entangled two-qubit states, re- 
quires (at least) one of the degenerate eigenvalues 1/4 of 
the partial transpose to become negative. This cannot 
happen within the perturbative regime. This behaviour 
can be readily understood in terms of the monogamy of 
entanglement, see, e.g., Ref. [24]. To first order in the 
small parameter expansion the reduced states of the ini- 
tially maximally entangled modes, A and k or n' and C, 
respectively, remain unperturbed up to relative phases 
due to the time evolution. In particular, the bipartite en- 
tanglement between A and well as between n' and C 
is maximal to first order in h, see Ref. [7], which excludes 
the possibility of first order correlations between k and 
k'. 

Furthermore, we can straightforwardly classify the en- 
tanglement of the four qubit state. We notice that, to 
first order, the state Q^kh'C can ^ e decomposed as 

QA KK ' C = mm\ + o(h 2 ) (l?) 

where || T> )) is a Dicke state [25, 26], given by 
|X>) = | ( | )) ± G* K |] Uf 1 1„ )) + ± G K i I l K , 1 1 laf 
+ G K ,AfJ ||U + \1 K ,)T - GIAZ ||U)n|lc» + 

+ G:G K -iu)nL» + iw)riic» + ), m 

where G n are mode dependent phase factors of unit 
magnitude, i.e., \G n \ = 1, that arc determined by the 
specific travel scenario, see, e.g., Ref. [7]. The usual form 
of the Dicke state can be obtained from Eq. (17) by local 
unitaries, e.g., bit-flips in the modes A and k! ' . 



III. SINGLE-CAVITY ENTANGLEMENT 

Let us now modify our previous setup and consider 
only Rob's cavity on its own, as in Ref. [8]. In particular, 
let us assume for simplicity that the initial state of 
the in-region modes in that cavity is the vacuum. We 
want to investigate the genuine multipartite correlations 
that might possibly be generated between three selected 
modes by performing a Bogoliubov transformation to 
the out-region modes. 



A. Scalar fields 

The in-region vacuum | ) of a real scalar field <f> can 
be related to the out-region vacuum | ) by | ) = 
Me w | ) , where M is a normalization constant and 
W := \ jVija\a*y The coefficients Vij form a sym- 
metric matrix that can be expressed as V — — /3*a , 
where (3 = (/? mn ) and it is implicitly assumed that 
d = (a mn ) is invertible. We then apply the small pa- 
rameter expansion for the Bogoliubov coefficients and 
find V = V {1) + 0(h 2 ) and the normalization constant 
M = (1 - ± \Vij\ 2 ) + 0(h 3 ), see Ref. [8]. 

Having transformed the initial vacuum state | ) to 
the out-region, we continue by tracing over all out-region 
modes except three chosen modes k,k' and k" . We de- 
note the transformed, reduced state by pkk'k", where we 
keep terms up to second order in h. At this stage we 
further assume that the modes do not all have the same 
parity, e.g., let us choose (k — k') and (k' — k") to be 
odd, which implies (fc — k") is even. This further means 
that | 0£l | > 0, \0$ k „\ > 0, while \P k %,\ = 0. With 
this convention in mind we select a witness for genuine 
multipartite entanglement, i.e., 

2\{Q\p kk , k n |l fc )|2fc')|lfc")l - O(h 3 )<0, (19) 

which can be rewritten as 

2V2|/3< 1 »,||/3^„| - O(h*)<0. (20) 

As previously, the witness (19) presents a lower bound to 
the convex roof extension of the minimal average mutual 
information across all bi-partitions [21]. It can be imme- 
diately noticed that that the previously discussed bosonic 
resonances [10] allow the linear enhancement of the indi- 
vidual coefficients, 1/3^/1 or | |, for particular basic 
travel times r„ = 2nS/(k + k') and r m = 2mS/(k' + fc"), 
n,m <G IN+, respectively. Interestingly, these resonances 
coincide for n = p(k + k') and n = p(k' + k"), such that 
the travel time, as measured at the centre of the cavity, 
for a single cycle of the repeated basic travel scenario is 
t = 2pS, p G 1N+, see Fig. 2. 

At this mode-independent resonance the lower bound 
on the genuine multipartite entanglement increases 
quadratically with the number N of repetitions of the 
basic travel block. At the same time the mixedness that 
is introduced to the system by tracing out the other 
modes contains second order terms /f^ fe /, f^ kk „ and 

/f„^„ where f^ p = ±E„ #P «| 2 , which all ex- 
hibit quadratic scaling at the mode- independent reso- 
nance. The validity of the perturbative approach is en- 
sured since all second order terms are at most propor- 
tional to N 2 h 2 < Nh. 

However, the classification of the bosonic genuine 
multipartite correlations remains an unsolved problem. 
To first order the transformed state can be written as a 
pure Dicke state [26], but since the genuine multipartite 
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entanglement is detected at second order this is of no 
significance. To second order the modes effectively 
become qutrits, for which generally little is known about 
entanglement classes. 



B. Dirac fields 

Let us again consider the fermionic counterpart to the 
bosonic situation. The in-region vacuum || )) of the 
Dirac field tjj is related to the out-region vacuum || )) by 
I )) = Me w | 6 )), where W := £*>° V M &t c t an d M is 

a normalization constant. Working pcrturbatively in the 
parameter h the coefficient matrix V can be expanded in 
a Maclaurin series as V = V [1) +0(h 2 ) and V$ = G q A^ q * 
(no summation), where G n is a mode specific phase fac- 
tor, i.e., \G n \ = 1, that depends on the chosen travel 
scenario, see Ref. [7]. The normalization constant can be 
found to be M = 1 - \ >o |V$ | 2 + 0{h 3 ). 

q < 

We can then perform the Bogoliubov transformation 
of the in-region vacuum || )) to the out-region state and 
trace over all modes except three chosen modes K > 0, 
k' > and k" < 0, obtaining the state g KK ' K ". Since 
we do not expect any coherence effects between modes 
of the same charge when we start from the vacuum, see 
Ref. [8], we further assume that (k + k') is even, while 
(k+k") and (k'+k") are odd. Specializing to the massless 
case this implies that |*A^„| > and \As; K „\ > 0, while 

kd = o. 

For the state q kk ' k " wc can form a witness for gen- 
uine multipartite entanglement using the techniques from 
Ref. [27]. The violation of the inequality 

« I Q K H , K „ || 1 K »+ | 1„» ))- | + | « I Q K K , K „ 1 1 K , »+ | 1„» ))- | 

- [« o 1 g KK , K „ || o )) (iv \+(i n 1 l i n f | i H « )>- 

+1V || + «1 K < II Q KK ' K " II 1 K > t II L- >r)] * - ©(ft 2 ) < (21) 

detects genuine multipartite entanglement. It can be 
expressed in the simple form 

I- 4 ™" I + I - \/|^Lk''I 2 + I^«'k„| 2 + 0{h 2 ) < . (22) 

Using the triangle inequality this can be easily seen to 
be violated whenever both Ailn and ^4.„7„« arc nonzero. 

Wc can further classify the genuine multipartite entan- 
glement in this case, since the state Qk k,' k" admits the 
decomposition 

Qk«k* = \\W))(W\\ + 0(h 2 ) (23) 

where the class-defining W-state is 

||W))= ||6)) + G K „A%, \l K ')f\\l K '-f 

+ G K »A<%, IlLflVF • (24) 



IV. CONCLUSIONS 

We have studied genuine multipartite entanglement 
of bosonic and fermionic modes of rclativistic quantum 
fields in non-uniformly moving cavities. We have used the 
perturbative approach of Refs. [G, 7] to handle the Bogoli- 
ubov transformations that feature in the transformation 
of the cavity modes between the inertial in-region and 
out-region. The non-uniform motion in between these 
regions populates modes and shifts pre-existing excita- 
tions. The final out-rcgion states of a set of chosen modes 
is then obtained by tracing over all other modes. In two 
qualitatively different scenarios we have shown that gen- 
uine multipartite correlations are generated from initially 
bi-separable or separable states of the chosen modes. 

We have employed witnesses for multipartite entangle- 
ment that prove to be advantageous with respect to usual 
entropic measures in the perturbative regime. The wit- 
nesses for the bosonic correlations provide lower bounds 
to measures of genuine multipartite entanglement that 
are based on convex roof extensions of the minimal av- 
erage mutual information over all bi-partitions [21]. We 
find that the numerical value of the perturbative correc- 
tions to these lower bounds can be resonantly enhanced 
for any chosen triple of bosonic modes with varying par- 
ities. However, the classification of the arising correla- 
tions is hindered by the unknown classification structure 
beyond qubits. 

For the fermionic systems, on the other hand, wc 
can detect the genuine multipartite entanglement in the 
transformed states and, subsequently, assign the entan- 
gled states originating from initially bi-separable and sep- 
arable states to the classes of 4-qubit Dicke states, which 
can be used in quantum secret sharing [28], and 3-qubit 
lU-states respectively. 

The creation of specific entangled states in our setup 
can be considered as the realization of quantum gates by 
motion in spacetime. In particular, we complement the 
findings of Refs. [9, 10], where two-mode squeezing gates 
are implemented as a result of the non-uniform cavity 
motion, with the creation of the Dicke and W states for 
fermionic systems. 

Indeed, the cavity setups studied here and in Refs. [8- 
10] share intriguing features with the models for fre- 
quency combs, which are known to produce cluster states, 
a vital resource for universal quantum computation [29]. 
Future work is being directed towards the investigation of 
this connection as well as to the extraction of the cavity 
mode entanglement with suitable detector models [30]. 

Additionally, the presence of the genuine multipartite 
correlations in these relativistic settings can be used for 
high precision tests of the quantumness of correlations 
and might have significant advantages in identifying the 
signatures of quantum phenomena where relativistic ef- 
fects are notoriously small. The resonant behaviour of 
the bosonic multipartite entanglement presented here can 
be a keystone countermeasure to this problem and al- 
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lows precise control and, hopefully, experimental testing 
of such effects. 
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